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TtO08 . . IST eS See TTT ATTA TC a 
-TITLE MTHSGTANH : A jos ting Hyperbolic Tangent routine 
IDENT /1-005/ + File: MTHGTANH.MAR Edit: JCw1005 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
L EQUIPHENT CORPORATION, MAYNARD, MASSACHUSETTS. 
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afACdLitys MATH LIBRARY 
ABSTRACT: 
MTHSGTANH is a function which returns the G pet hte hyperbolic tangent 


of its G floating point argument. The call is standar 
tall- by-reference. 
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; Detailed Current Edit History 6=SEP-1984 MTHRTL.SRC at THGTANH.MAR; 1 (2) 


-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
Edit History for Version 1 of MTHSGTANH 


1-001 = Adapted from MTHSDTANH 24esepter3 1- 7906. SBL 26-Jan-79 
- ¢ Use MTH$GEXP_R6. opt Be 
1-005 - enone constant 16 © 90, ty correct inaccuracy. The 
e of X above ae is the best machine nbers thag cen 
to GTANH(X) is about 19. 08. The oe | sina ous r that can 
be represented as a short literal -Sep-80 
Eliminated symbolic short gtisecals “Suh j coct-81 
Changed the constant 2“-1 7 to correct inaccuracy. A 
values of 1X! between pa.f2 and 2*=-27 GTANH was only accura 
he 8 decimal preges since the gesagt ion that hag bday Ay that 
onge of values is false. ap aay for iXi<=2*=27 a 
GTA H(X)= GSINH(X) /GCOSH(X) for 2%-c7<ixi<=.25 . ALL appropriate 
references to 2“-14 have been changed to 24-27, JCwW Jan-83 
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CLARATIONS DB g: E 71982 ot: 4:27 ° | 


eclarative Part of Modul SEP MTHRTL.SRCJMTHGTANH.MAR; 1 
»SBTTL DECLARATIONS ; Declarative Part of Module 


: INCLUDE FILES: 
; EXTERNAL SYMBOLS: MTHSJACKET_HDLR 


-DSABL GBL ; Force .EXTRN on all symbols 
~EXTRN MTHSGCOSH ; GCOSH 
»EXTRN MTHSGSINH ; GSINH 
~EXTRN MTHSGEXP_R6 3 EXP 
; EQUATED SYMBOLS: 
value = 4 ; value.rg.r 
> MACROS: none 
> PSECT DECLARATIONS: 
-PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 


3; program section for math routines 
OWN STORAGE: none 


CONSTANTS: 
G_0.25: 
-WORD “X3FFO, 0, 0, 0 ; 0.25 
G_2_POWER M27: 
-QORD “X3E60, 0, 0, 0 3 2e8=27 


a 
MEW SO OO NA UE WN O ODNOA UE WWN COONS Ur 
. 


SGOOCOCCOCOCOOOOOOOOODWOWOOWDOMWMOO IV Y~ 


¢ 9 


MTHSGTANH G Floating Hyperbolic Tangent routine 16-SEP-1984 01:32:44 VAX/VMS Macro v04-00 Page 4 
1-008 MTHSGTANH : G Bouble Precision Floating att f 0:32:58 EMTHRTL. SRE IMTMGTANH MARS 1 ° (4) 
! 108 -SBTTL MTHSGTANH = G Double Precision Floating DTANH 
10 109 
1 110 ;++4 
! 1} ; FUNCTIONAL DESCRIPTION: 
a ae ; GTANH = G floating point function 
: 12 ; GTANH(X) is computed as: 
61 119 ; If iki =< 2**-27, then GTANH(X) = X. 
Bat 113 3 If 2e*-27 < {Xi = -25, then GTANH(X) = GSINH(X)/GCOSH(X). 
10 119; If 0.25 < !X! < 20.0, then GTANH(X) = (GEXP(2*x) = 1) / (GEXP(2*x) + 1) 
3018 0 ; If 20.0 =< !X!, then GTANH(X) = sign(X) * 1 
001 1 § 3; CALLING SEQUENCE: 
010 123; 
0010 124 ; GTANH.wg.v = MTHSGTANH(x.rg.r) 
0010 125; 
Bo10 : 8 3 INPUT PARAMETERS: 
00000004 0010 128 LONG = 4 ; define Longuorg multiplier 
00000004 Baio : ; x = 1 * LONG ; Contents of x is the argument 
0010 131 ; IMPLICIT INPUTS: none 
0010 13 ; 
0010 133 ; OUTPUT PARAMETERS: 
0010 134 ; 
O38 133 ; VALUE: G floating hyperbolic tangent of the argument 
0010 it ; IMPLICIT OUTPUTS: none 
0010 138; 
0010 139 ; COMPLETION CODES: none 
0010 140 ; 
0010 141 ; SIDE EFFECTS: none 
0010 146 3 
0010 1435 ; NOTE: This procedure disables floating point underflow, enables integer 
0010 144 ; overflow. 
0010 145 ; 
0010 146 gore 
0010 14 
001 148 
407C 001 149 -ENTRY MTHSGTANH, “M<IV, R2, R3, R4, RS, ROD> 
oot 150 ;standard call-by-reference entry 
01 151 ; disable DV (and FU), enable IV ; 
Bat 152 MTHSFLAG_JACKET : flag that this is a jacket procedure in 
6D Q0000000'GF 9E ot MOVAB G*MTHSSJACKET_HND, (FP) 
01 ; set handler address to jacket 
B12 3; handler 
619 153 3; case of an error in routine 
19 154 ; If an error, convert signal to user PC 
19 155 : and resignal 
50 gs C SOFD 19 136 MOVG Oya lye tAr? RO ; RO/R1 ‘Xi = @value(AP) 
50 8000 8F AA Ole 1 Bicw  #*x8000, RO : RO/RI = 1X} 
EO AF 50 51FD 0023 158 CMPG = RO, -G_2_POWER_M27 : compare !X! with 2**-27 
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BLEQ 
2ee=27 < 1X! 


CMPG 
BGEQ 


ouT_X 


RO, #20 
GE4_T0_20.0 


; 2ee=27 < 1K! < 20.0 


CMPG 
BLEQ 


RO, G.0.25 
LE4_10_0.25 


; 0.25 < ixt < 20.0 


ADDG3 


; 2ee-1R6 < IX! 
LEQ_T0_0.25: 
CALLG 
V 
ALL 
Vv 


CALLG 
DIVG2 
RET 
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GEQ_T0_20.0: 
MOVG 


TSTG 
BGEQ 
MNEGG 
RET 


@value(AP), a@value(AP), RO 
W*MTHSGEXP_R6 
1 R2 


R2, RO 


=< 0.25 


(AP), W°MTHSGCOSH 
(AP), W*MTHSGSINH 
R2, RO 


#1, RO 
@value(AP) 
10$ 

RO, RO 


@value(AP), RO 
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; compare 
; branch if 


; compare 
; branch if 
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; RO/R1 = GCOSH(X) 
; R2/R3 = 
; RO/R1 = GSINH(X) 
GSINH(X) / GCOSH(X) 
return with result in RO 


RO/R1 = 


RO/R1 = 1.0 
test the sign of x 


branch if 
RO/R1 = 


return with result in RO 


RO/R1 = GTANH(X) = X 
return with result in RO/R1 
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Symbol table ~ petite Epayo8e S1i2ese? EMTRTE eae oaTNGtANc.maR:1 2% chy 
Gea -J 20.0 0 R 1 | 
g-9¢ R 1 
G- R_M27 R 1 
LEo~ 10, 25 R 1 
LON 0000004 
RINSSJACKET _HND teeneresr x 4 
MTHSGCCOSH eereeeee x 0 
MTHSGEXP_R6 eeeeeree 44 
MTHSGSINA eereeene x 
MTHSGTANH 00000010 RG 4 
OuT x 00000071 R 1 
VALOE = 00000004 
peewee aocsneeeewe roe + 
! ; _Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes | 
- 4ABS. 00000000 0.) 00 ¢ OO.) NOPIC USR CON ABS LCL _*- NOEXE NORD NOWRT NOVEC BYTE | 
_MTHSCODE 00000076 (¢ 118.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG | 
foewoccoeceecesooeceseccena} 
H Performance indicators : 
Phase Page faults CPU Time Elapsed Time | 
Initialization | 30 0880 00-08 69:00:00.7% 
zonnene processing | OF B88 82 39:80:89 -02 
ss :00:00. 00:04. 
pyases table sort 59 00:00: 3°9 00:00:00-00 
ass H : e : H . 
Symbol table output ; 00:00: §°3 $000,008 
Psect synopsis output 00:00:00.0 bs SS 
Cross-reference output ° 00:00:00.0 8:99:09 00 
Assembler run totals 28 00:00:02.06 00:00:10.31 


The working get Limit was 900 pages. : 

3309 bytes (7 pages) of virtual memory were used to buffer the intermediate code. 

aed were 10 pages of symbol table space allocated to hold 13 non-local and 1 local symbols. 
278 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to “define 1 macro. 


weecew ee ses ese ee eo sa eee eee = + 


Macro Library name Macros defined 

“S255SDUAZB:(SYSLIBISTARLET.MLB:2 0 tS 0 

0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. | 


MACRO/ENABLE=SUPPRESS1ION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=L1S$:MTHGTANH/OBJ=OBJ$:MTHGTANH MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRC | 
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